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Abstract 


This paper is devoted to the derivation of simple approximate equations for transmission and reflection 
functions of multi-layered disperse media. It is assumed that each layer is optically thick and weakly 
absorbing. Equations obtained can be used in various fields of applied optics, including cloud and snow 
optics. 
© 2004 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Vertically inhomogeneous disperse media are of a frequent occurrence both in nature (e.g., 
multi-level cloud systems, snow deposited at different times at a given place, terrestrial 
atmosphere and ocean, biological tissues, etc.) and technological applications (multi-layered 
painted surfaces, paper, etc.). This explains a great interest in studies of radiative transfer in 
vertically inhomogeneous media. Recent advances in this area have been summarized by 
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Yanovitski [1]. A great number of exact (see, e.g., [2,3]) and approximate [2,4-7] techniques has 
been developed. At present time there is no problem to account for an arbitrary vertical 
inhomogeneity of a horizontally homogeneous plane-parallel light scattering medium. However, it 
appears that for this, one needs to perform quite complex numerical calculations with the use of 
up-to-date computer technology. On the other hand, the practical work requires simple 
approximate solutions, which can be used to perform rapid estimations of the influence of vertical 
inhomogeneity on light reflection and transmission by turbid media. With this in mind, we have 
decided to derive such approximations based on our recent work for a single homogeneous 
optically thick disperse layer [8] and also on earlier works of Rozenberg [9], Melnikova and Minin 
[10], and Zege et al. [11]. To derive approximate solutions, we make an assumption that the 
probability of light absorption by particles in a disperse medium is close to zero (smaller than 
approximately 3%). 


2. Theory 
2.1. Reflection 


Let us assume that a disperse medium is composed of several light scattering layers having 
different light scattering characteristics (e.g., the phase function p(0), the absorption Gabs and 
extinction dex; coefficients). The geometry of the problem is shown in Fig. 1. Light enters a 
disperse medium at the angle 0o. The reflected light is observed at the angle 0; and the diffusely 
transmitted light is observed in the direction specified by the angle 02. We need to model the 
intensity of diffusely reflected and transmitted beams taking into account effects of vertical 
inhomogeneity of the medium under consideration. 

Note that we put no limitations on the number of layers in a disperse medium. In Fig. 1 we 
show a situation, which can be a case in the terrestrial atmosphere with several cloud layers at 
different heights. Then the height in Fig. 1 is given in kilometers. So we refer to these layers as 
clouds in this paper. Although layers can be composed of other disperse substances and materials, 
the region which separates clouds is assumed to be free of light scattering and absorption. In 
particular, it means that the situation we consider is equivalent to the touching light scattering 
layers (e.g., several snow layers). This assumption, however, is not crucial for the theory 
developed here (e.g., absorption and scattering of light by gases and aerosols between clouds can 
be accounted for, if needed (see, e.g. [8])). We assume that all layers are optically thick (L;>/;, 
where L; is the geometrical thickness of i-layer and l; is the photon free path length in the layer, 
which is equal to the inverse value of the extinction coefficient in a given layer) and weakly 
absorbing (the probability of photon absorption p = Gabs/Gext > 0). We need these assumptions 
to apply the modified asymptotic theory for each layer as specified by Kokhanovsky and Rozanov 
[8] (see Appendix A). There are no limitations on the type of the phase function, however. 
Although our assumptions severely restrict the applicability of the model to many real-world 
media, they do provide an accurate approach to the solution of a number of important problems 
(e.g., in optics of clouds, snow, biological tissues and paints). 

The starting point is the expression for the reflection function of a single homogeneous optically 
thick weakly absorbing layer. It can be written in the following form for a scattering layer above a 
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Fig. 1. The geometry of the problem. 


Lambertian surface with the albedo A [12]: 
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1 — Ar () 


Ralé, 1, 1) = RE, N, 1) + 
Here (see Appendix A) 
RE, n, T) = RxlE, n) z Ie, 1,7) exp(—x =y) (2) 


is the reflection function of a scattering layer for the black (A = 0) underlying surface, R4(é, n, T) is 
the same function but for the scattering layer—underlying surface system, x = kt, k = \/3fB(1 — g) 
is the diffusion exponent, g is the asymmetry parameter, which is easily obtained by the 
integration of the phase function with respect to scattering angle [13], tT = Gextl is the optical 
thickness, / is the geometrical thickness of a scattering layer, y = 4k/3(1 — g), r is the spherical 
albedo of a scattering layer for an illumination from below at A = 0, ¢g(€) is the diffused 
transmittance for the illumination in the direction 09 = cos“! (é) (see Fig. 1), T(é,7,7) is the 
transmission function at A = 0. The function ta(n) is the diffused transmittance for the 
illumination in the direction 0, = cos~!(y). The dependence of the reflection function in Eq. (1) on 
the azimuth ¢ is omitted for simplicity. Also we neglect the direct light transmittance, which takes 
rather small values for optically thick layers. Convenient approximate equations for functions 
T(é,,7), ta(E), Roo(€,4) and r are given in Appendix A (see also ref. [13]). 
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Let us proceed further now. It is known that the reflection of light from an optically thick 
weakly absorbing strongly light scattering layer is rather close to that of a Lambertian reflector. 
This is due to well-developed multiple light scattering in the medium in this case. Therefore, to 
find the reflectance from a multi-layered system as given in Fig. 1, we propose to use Eq. (1) with 
A substituted by the spherical albedo r= of the multi-layered system below the upper layer (see 
Fig. 1). Other functions in Eq. (1) then refer to the first layer from the illumination side. Such an 
approach was proposed and successfully used by Melnikova and Minin [10] for studies of light 
fluxes in a cloudy atmosphere. 

The spherical albedo r5 can be easily found integrating Eq. (1). Indeed, it follows from Eq. (1) 
(see Appendix A): 


= (3) 


Here r¥ is the spherical albedo of the system starting from the third layer down (see Fig. 1), f is 
the global transmittance of the second layer (see Appendix A for the definition), rz is the spherical 
albedo of the second scattering layer for an illumination from below. The spherical albedo r¥ can 
be found using an equation similar to Eq. (3) 


z (4) 

with meaning of all parameters similar to those in Eq. (3). Clearly, we need to repeat this 

procedure till the underlying surface is reached. Then we have 
A 

1—r,A 


r =n + (5) 
and the procedure is complete. 

Let us check the accuracy of the straightforward procedure outlined above using the exact 
solution of the radiative transfer equation with SCIATRAN [3] for a two-layered disperse system 
over a black surface. Then A in Eq. (1) should be substituted by the spherical albedo of a lower 
layer r2. All other parameters in Eq. (1) refer to an upper layer. Note that SCIATRAN is a well- 
documented and thoroughly tested radiative transfer code based on the discrete ordinates 
approach. Its accuracy is better than 1%. 

The results of comparisons are shown in Figs. 2a—c. In particular, we give the dependence of the 
reflection function of a disperse medium on the incidence angle for the nadir observation. Due to 
the reciprocity principle, our calculations are also valid for the nadir illumination and varying 
observation angles. 

The middle curves in Figs. 2a—c correspond to a two-layered cloud system with the optical 
thickness of the bottom layer tp = 30 and the optical thickness of the upper layer ty = 10. We 
assume that the upper layer does not absorb incident radiation. The single scattering albedo of the 
bottom layer is equal to 0.9945 (Fig. 2a), 0.9982 (Fig. 2b) and 0.9681 (Fig. 2c). 

The upper curves in the same figure correspond to a single layer having optical thickness 
T = Tb + Tu and the single scattering albedo equal to 1. Lower curves in Figs. 2a—c correspond to a 
single layer having optical thickness t = Tp + Tu and the single scattering albedo equal to 0.9945 
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Fig. 2. (a) The dependence of the reflection function of a cloud medium on the incidence angle for the nadir 
observation. Symbols give exact results and lines are due to the approximation (see details in text). The curve in the 
middle was obtained for an upper nonabsorbing layer and the absorbing layer with wo = 0.9945 at the bottom. Upper 
curves correspond to a nonabsorbing single layer. Lower curves correspond to an absorbing single layer. The total 
optical thickness is kept constant for all calculations (t = 40). (b) The same as in (a) except wọ = 0.9892. (c) The same as 
in (a) except wp = 0.9691. 


(Fig. 2a), 0.9982 (Fig. 2b) and 0.9681 (Fig. 2c). Exact results are shown by symbols. Lines 
correspond to calculations according to the approximation developed here. 

The phase function was found using the narrow gamma droplet size distribution [13] with the 
effective radius 10 um and the effective variance equal to § at the wavelength 2 = 0.65 um for 
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water droplets with the refractive index m equal to 1.331 — 0i (Figs. 2a—c, upper curves), 1.331 — 
0.00005; (Fig. 2a, lower curves), 1.331 — 0.00017 (Fig. 2b, lower curves), 1.331 — 0.0003: (Fig. 2c, 
lower curves). 

Note that the variation of the imaginary part of the refractive index allows us to model various 
levels of cloud pollution (e.g., due to black carbon). The phase function differs not significantly 
for all layers. For instance, the asymmetry parameter equal to 0.85 at m = 1.331 — 0i, 0.8513 at 
m = 1.331 — 0.00005i, 0.8525 at m = 1.331 — 0.00017, and 0.8571 at m = 1.331 — 0.00037. This is 
according to the general fact that the phase function of weakly absorbing particles is not 
particularly effected by the level of light absorption. However, note that there is a slight tendency 
for the general increase of the asymmetry parameter with the imaginary part of the refractive 
index. 

It follows from Figs. 2a—c and 3 that the accuracy of our simple approximation is quite high for 
the case considered. In fact it is comparable with the accuracy of correspondent equations for a 
single layer or even higher than that (e.g., see Fig. 3 for a single layer at wo = 0.9691). This 
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Fig. 3. The relative errors 6 = 1 — R4/Re (in percent) for calculations shown in Fig. 2. Here indices a and e signify 
approximate and exact results, respectively. 
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paradox is explained by the fact that the accuracy of the approximation is highly influenced by the 
value of wo. The average value of wo is, however, lower for a two-layered system (with a 
nonabsorbing upper layer) as compared to a single absorbing layer having the same optical 
thickness. 

We see that the two-layered system with the total optical thickness 40 has values of R 
intermediate between those for an upper layer (at t = 40) and lower layer (also at t = 40). This 
can be expected on general grounds as well. Errors are generally below 5% but they increase for 
oblique incidence angles. The accuracy decreases with f. The value of f ~ 0.03 can be considered 
as an upper boundary for the application of this theory. Although it can be applied to slightly 
larger values of f if the accuracy is not a primary concern (e.g., for rapid estimations of vertical 
inhomogeneity effects). 

It is interesting to see the performance of equations for larger and smaller values of tp, ty. This 
is shown in Figs. 4 and 5 for an absorbing lower layer and nonabsorbing upper layer. We see that 
the accuracy of our equations is robust against change of the turbid layer thickness. Note that the 
variation of the optical thickness of a lower layer (see Fig. 4a) does not change the reflection 
function very much. This is due to the fact that the spherical albedo of a lower absorbing cloud 
does not depend strongly on ty. On the other hand, the variation of the upper layer optical 
thickness Tu (see Fig. 5a) changes the result considerably. Cloud becomes much brighter for the 
larger thickness of the upper layer. Obviously, for a very thick upper layer the sensitivity of the 
reflection function to the presence of a turbid layer at the bottom is lost. This is similar to 
the effect of the disappearance of objects in a heavy fog. Therefore, we conclude that high 
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Fig. 4. (a) The dependence of the reflection function of a two-layered cloud medium on the incidence angle for the nadir 
observation at Tu = 10 and tp = 5, 10, 30. The upper layer does not absorb radiation. The bottom layer is characterized 
by the single scattering albedo 0.9892 (see Fig. 1b). Symbols give exact results and lines are due to the approximation 
(see details in text). (b) The relative errors 6 = 1 — R4/Re for calculations shown in (a). 
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Fig. 5. (a) The dependence of the reflection function of a two-layered cloud medium on the incidence angle for the nadir 
observation at tT, = 5, 10, 30 and tẹ = 30. The upper layer does not absorb radiation. The bottom layer is characterized 
by the single scattering albedo 0.9892 (see Fig. 1b). Symbols give exact results and lines are due to the approximation 
(see details in text). (b) The relative errors 6 = 1 — R,/R- for calculations shown in (a). 


nonabsorbing clouds can shield lower (and possibly) polluted clouds. This can lead to important 
climatic effects not accounted in Global Circulation Models at the moment. 

To make this point more clear we present the reflection function of a single absorbing cloud 
layer with the optical thickness 30 in Fig. 6. Then we add a nonabsorbing cloud at the higher level 
in the atmosphere. It follows that the reflection of the system considerably increases both for 
warm water and cold ice upper-level clouds. The phase function of an ice cloud was taken from 
study of Mishchenko et al. [14] (the fractal particle model) and the phase function of the water 
cloud was calculated as indicated above (at m = 1.331 — 0.00017 for a lower cloud and m = 1.331 
for an upper cloud in the two-layered system). The increase in the reflection is much more 
pronounced for crystalline clouds. It means that ice clouds not only warm the system by trapping 
terrestrial radiation. They also may shield lower polluted cloud systems (e.g., in urban areas) and 
increase general reflection of the surface—atmosphere system. This indicates the complexity of the 
issue of cloud’s influence on climate. 


2.2. Transmission 


The question arises if the model presented above can be applied to studies of light transmission. 
The answer to this question is positive. Indeed, we have for the transmission function of a single 
homogeneous disperse layer over a Lambertian surface [12] 


Atal RC) 


Tal, n, T) = T(&, n, 7) + FT 


(6) 
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Fig. 6. The dependence of the reflection function on the incidence angle for the nadir observation for a two-layered 
cloud medium at Tu = 10, tp = 30 (upper curves) and for an absorbing single cloud layer with t = 30 (lower curves). 
The upper layer does not absorb radiation and is composed either of water droplets or ice crystals. The bottom layer 
and a single absorbing layer is characterized by the single scattering albedo 0.9892 (see Fig. 1b). Symbols give exact 
results and lines are due to the approximation (see details in text). 


Here all functions are defined as in Eqs. (1) and (2) and in addition R(n) is the plane albedo of the 
layer for the illumination from below (see Appendix A). Again assuming that the layers below the 
upper one can be substituted by a Lambertian reflector, we have for the transmission function 
T\(é,n) of the first layer in the n-layer system 


ta(E)Ri(n)r5 


T\(é.n) = T1(, n) + T= riž 
2 


i (7) 
where we omitted the dependence on the optical thickness and rž is found using the iterative 
procedure starting from the ground surface as underlined in the previous section. Functions 
T\(é,7), tai(€), rı and (4) have the same meaning as in Eq. (6) but for the first layer. 

Let us consider now the transmittance under the second layer. The second layer is illuminated 
by the diffused light transmitted by the first layer. It follows from Eq. (7) that the diffused 
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transmittance 74;(€) (see Appendix A) is given by 


t abe 
uO = u + 8 (8) 
= rir} 
or 
-oan talë) 
WO = ape ©) 
Also we have for the global transmittance from Eq. (9) (see Appendix A) 
a ti 
Îi DAF <7 (10) 


Therefore, the transmission function T, after the second layer in the n-layered system is given as 


T(E.) = ta (€)ta2(), (11) 
where (see Eq. (9)) 
fat) = po. (12) 
3 


Following this procedure, we can obtain the transmission function under the third layer 


T3(E,n) = fa (&)h ta(n), (13) 


where 7) = t2/(1 — r2r3) and we accounted for the fact that the second layer is illuminated from 
above by diffused light and it also serves as a diffuse light source for the third layer. Repeating this 
procedure for each layer, we can arrive finally to the transmission function of a whole system 


T(E.) = tai (Qiks... in_1Fan(), (14) 
where 7 = ¢;/(1 — rey) and r,,; = A. Interestingly, Eq. (14) can be written in the form similar to 
that for a homogeneous layer (see Appendix A) 


T(E,n) = ter Ko(€)Ko(n), (15) 
where the effective global transmittance is given by 
I- 1G 


(16) 


lef = a 
Ip- 0 — ya) 


with all parameters defined in the theory for a single layer. Note that we have used here the 
equality: tg; = t;Ko(¢), where Ko(¢) is the escape function (see Appendix A). 

Let us check the applicability of our assumptions making comparisons with exact radiative 
transfer calculations using SCIATRAN [3] for a special case of a two-layered medium over a 
black surface. Then Eq. (15) reduces to the following form: 


tt2Ko(¢)Ko(n) 
1— rira 


T(E, n) = (17) 


where we accounted for the fact that tet = t)f(1 — rjr2)7! in this case. Note that if both layers do 
not absorb radiation, then the sensitivity of transmitted light to the vertical inhomogeneity is low 
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and in a good approximation one can use the reflection function for a single layer having the 
optical thickness equal to the sum of optical thicknesses of both layers and the average value of 
the asymmetry parameter to find the solution of the problem at hand [12]. 

The results of calculations using simple approximation (17) are shown in Figs. 7—9 as functions 
of the observation angle together with correspondent errors and outcome of exact calculations for 
the incident angle equal to 60° and the azimuth equal to 0°. In particular, we give the dependence 
of the transmission function on the observation angle for a single layer having the optical equal to 
40 at wo = 0.9945 and 1.0 in Fig. 7a. The results of computations for a two-layered medium 
having the total optical thickness 40 but mp = 0.9945 in the bottom layer (tp = 30) and wp = 1.0 
in upper layer (ty = 10) are also given in the same figure (the middle line). Note that phase 
functions in all calculations given here are very close to each other as specified above. So the 
change of transmission is mostly due to the absorption effect. As expected the largest transmission 
is observed for a nonabsorbing single layer. It is considerably reduced if absorption is introduced 
in the bottom part of a layer. Of course, the minimum of transmission occurs for a single 
absorbing layer (see, e.g. a lower line in Fig. 7a). It follows that exact and approximate results are 
quite close to each other for observation angles smaller than 70°. Then the error of approximation 
is smaller than 5% (see Fig. 7c). The error increases for slabs having larger absorption, however. 
This is illustrated in Fig. 7b, where we show the dependence similar to that in Fig. 7a but for the 
increased absorption (wọ = 0.9892). It is interesting that the error of Eq. (17) for a two-layered 
medium is smaller than that for a single layer with the optical thickness 40 and wp = 0.9892. This 
points to the fact that the accuracy of our technique is mostly influenced by the total light 
absorption and transmission and not by a number of layers (see also Fig. 7c). Note that the 
accuracy of reflected light calculation is generally higher than that for the diffusely transmitted 
light (for a given level of absorption). 

Calculations for a two-layered turbid slab with the optical thickness of a lower absorbing layer 
equal to 30 for various thicknesses of a nonabsorbing upper layer ta = 5, 10,30 are shown in 
Fig. 8a. The middle lines in this figure coincides with the middle line in Fig. 7a. We see that the 
accuracy is better than 5% in this case. The decrease of the optical thickness of an upper layer 
leads to the increase of the light transmission as one might expect. 

Fig. 8b is similar to 8a, but now the optical thickness of an absorbing layer at the bottom 
is varied from 5 to 30. The optical thickness of an upper layer is fixed and equal to 10. Clearly, 
the error approximation increases for thinner layers, which is in accordance with general 
assumptions of our approximation, which is valid only for weakly absorbing optically thick layers 
[15]. However Eq. (17) has comparatively high accuracy even at such small values of t as 5 (see 
Fig. 8c). 

In conclusion, we show the transmission function of a single absorbing layer having optical 
thickness 40 at wo = 0.9945 in Fig. 9a (middle line) in comparison with transmission functions of 
a two-layered system having wo = 0.9945 in the upper layer (ty = 30) and mp = 0.9892 in the layer 
at the bottom (tp = 10) (lower line in Fig. 9a). Clearly, the transmission for the latter case should 
be lower. Fig. 9a quantifies this decrease. Upper line in Fig. 9a corresponds to a two-layered 
system with total optical thickness equal to 40 and local optical characteristics of a lower layer 
equal to that of a single layer shown by the middle line in Fig. 9a but having a nonabsorbing 
scattering layer at the top of system (t, = 10). Then due to the general decrease of absorption 
in the system, transmission should increase. This is confirmed by Fig. 9a. It follows from Fig. 9b 
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Fig. 7. (a) The dependence of the transmission function on the zenith observation angle at the zenith incidence angle 
equal to 60° and the azimuth equal to 0° for a single homogeneous layer and a two-layered turbid medium with the total 
optical thickness equal to 40. The single scattering albedo is equal to 0.9945. Lines correspond to Eq. (17) and symbols 
are obtained from exact calculations. Further explanations are given in text. (b) The same as in (a) but at wo = 0.9892. 
(c) Errors of approximation (17) for cases considered in (a) and (b) (1, 2, 3—single layers with the optical thickness 40 
and wo = 1.0, 0.9945 and 0.9892, respectively; 4,5—two-layered systems with a nonabsorbing top layer having optical 
thickness 10 and an absorbing layer at bottom with ap = 0.9945 (4), wp = 0.9892 (5)). 


that the error of Eq. (17) is smaller than 5%, which is acceptable for a broad range of 
applications (e.g., rapid estimations of stratification effects on diffusely transmitted and reflected 
light fields). 
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Fig. 8. (a) The same as in Fig. 7a but for other values of optical thickness of an upper layer. (b) The same as in Fig. 7a 
but for other values of optical thickness of a bottom layer. (c) The error of approximation given by Eq. (17). 


3. Conclusion 


We derived here simple analytical equations for reflection and transmission functions of multi- 
layered turbid media in the assumption that each layer is optically thick (t>5) and weakly 
absorbing (f<0.03). They can be used to solve a variety of practical problems including light 
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Fig. 9. (a) The dependence of the transmission function on the zenith observation angle at the zenith incidence angle 
equal to 60° and the azimuth equal to 0° for a single homogeneous layer (middle line, wọ = 0.9945, t = 40) and a two- 
layered turbid medium (mp = 0.9945, t = 30 for a lower layer and wo = 1.0, t = 10 for an upper layer (upper line) and 
wg = 0.9892, t = 10 for a lower layer and mp = 0.9945, t = 30 for an upper layer (lower line)). Lines correspond to 
results obtained with Eq. (17) and symbols are obtained from exact calculations. Further explanations are given in text. 
(b) Errors of Eq. (17) correspondent to data given (a). 


propagation in stratified snow and clouds. The error of equations for the cases studied here is 
below 15% (for most of geometries the error is below 5%). 

Note that equations presented above should be used only in the case of absorbing layers. For 
nonabsorbing layers the sensitivity of reflection and transmission to the stratification is low. Then 
the theory of a single layer with an averaged asymmetry parameter and integral optical thickness 
equal to the sum of optical thicknesses of all layers should be preferably used. 
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Appendix A 


Here we present in a condensed form main equations for a single scattering layer proposed by 
Kokhanovsky and Rozanov [8]. Namely, it follows: 


RE, 1, 1) = RoE, n) = T(&, 1, 1) exp(—x E y), 
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where 


RoolEn, 6) = R? (E, n, p) exp(—(1 — 0.05y)u), 


sh(y) 


GE (=e + ay) 


=f; apa) Ko(&)Ko(n), 


a — bën + en E = Ko(€)Ko(n) 
T? ? R(E, p) 


KOS U+), te= 


_ ie 
x= V-D, y= Alaa, 


a = 4.86, b = 13.08, c = 12.76, f is the probability of photon absorption and g is the asymmetry 
parameter, č and 7 are cosines of incident and observation zenith angles, respectively, ¢ is the 
azimuth. The function R? (é, n, p) gives the reflection function for a semi-infinite nonabsorbing 
medium. It only weakly depends on the size of particles but considerably depends on their shape 
(see Fig. Al). We used the following parameterizations for this functions valid for the nadir 
observation geometry: 


0.37+1.94€ p(n — arccos €) 
1+¢ A(1 +é) 


for water clouds and 


R? (é, 1) = 


R? (E, 1) = 0.61284¢ + 0.8€ — 0.338492 


for ice clouds. The accuracy of these approximations is quite high as follows from Fig. A.1. 
Results in Fig. A.l are obtained for the gamma size distribution of water droplets f(a) = 
Ba’ exp(—9a/aer), B = const, a is the radius of droplets. The phase function of ice particles was 
taken from data given by Mishchenko et al. [14] for a fractal ice particle. The wavelength is equal 
to 0.65 um. 

The equation for the function R°(é,,¢) for other observation conditions is given by 
Kokhanovsky [16]. The phase function p(@) can be calculated using Mie theory. However, we used 
here the following approximate expression valid for water clouds [17]: 


5 
p(0) = ne™ + ` je 81-81)" 


i=] 


where 0 is the scattering angle in radians, n = 17.7, v = 3.9 and constants y;, 6;, 0; are presented in 
Table A.1. 

The reflection and transmission functions given above can be used to calculate various other 
radiative characteristics including the diffuse transmittance tą(¢), the global transmittance f, the 
plane albedo ‘R(€), and the spherical albedo r. Correspondent equations are given in Table A.2. 
Note that to derive formulae in Table A.2 from equations given above we neglected a small term 
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Fig. A.1. The dependence of reflection functions of a semi-infinite nonabsorbing water and ice cloud on the incidence 
angle at the nadir observation. Symbols give exact results and lines are due to the approximation (see details in text). 


Table A.1 

Parameters y;, 6;, and 0; [18] 

L Yi Ôi 0; 
1 1744.0 1200.0 0.0 
2 0.17 75.0 2.5 
3 0.30 4826.0 T 
4 0.20 50.0 T 
5 0.15 1.0 T 


te ~ t >. Also we used the equality 


1 
2 K dé = 1 
f (EE dé 


and approximate integration techniques described by Kokhanovsky [17]. 
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Table A.2 

Radiative characteristics of optically thick weakly absorbing strongly light scattering turbid layers [17] 

Radiative characteristic Symbol Definition Equation 

Plane albedo RE) 2 fÈ RE mndn Ral = Co — DKE) 
Spherical albedo r 2 H REE dé Too — t exp(—x — y) 
Diffused transmittance talé) i T(é, nn dn tKo(€) 

Global transmittance t 2f taé dé aeo 


Functions given in table are defined using following equations: 


R 


1 27 E E 1 2n 
af Ended T=} [Tene dd, 


3 
Roo(S) = exp[—W1 — 0.05y)Ko(O)], Kol) = 5 (1 + 26), 


Po =€, x=kt, yp=kt, k= V3(1—a@)l—g), «= 1.07. 


The light absorptance by the disperse medium under consideration for the mono-directional 
illumination A(é) and for the diffuse illumination a can be also derived. Namely, it follows using 
the energy conservation law: A(é) = 1 — ta(€) — R(é) and a = 1 -r — t. 
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